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Abstract 

We apply the factorization technique developed by Kuru et. al. [Ann. Phys. 323 (2008) 413] to obtain the exact 
analytical classical trajectories and momenta of the continuum states of the non Hermitian but VT symmetric 
Scarf II potential. In particular, we observe that the strange behaviour of the quantum version at the spectral 
singularity has an interesting classical analogue. 

Key words : Factorization method; Classical trajectories; VT symmetry; Scarf II potential; Spectral Singularity 

PACS Numbers : 03.65.-w, 03.65.Fd, 02.20.-a 

* e-mail : anjana23@rediffmail.com , sinha.anjana@gmail.com 



> 

00 
OV 

\6 
O 
(N 



X 



I. INTRODUCTION 

Non Hermitian Hamiltonians, especially those 
with VT symmetry or space time reflection sym- 
metry, aroused curiosity among the scientific com- 
munity ever since Bender et. al. showed that 
one can replace the condition of Dirac Hermitic- 
ity of the Hamiltonian (viz., = H) by the VT 
symmetry of the same (i.e. H VT = H), and yet 
obtain real and discrete spectrum under certain 
conditions [l], Q • Since then scientists have exten- 
sively studied such systems @, with an attempt 
to extend the framework of conventional quantum 
mechanics into the complex domain [4-11]. The- 
oretical predictions of such systems are found in 
quantum field theory, mathematical, atomic and 
solid state physics, classical optics, etc. 12j, while 



VT symmetric optical lattices have provided the 
ground for experimental verification [13-20]. 

The study of complex classical mechanics play 
an important role in understanding the classical 
limit of complex quantum theories, as it is gener- 
ally observed that quantum mechanics and classi- 
cal mechanics provide profoundly different descrip- 
tions of the physical world in many cases. In an 
attempt to bridge this gap, several studies (both 
analytical and numerical) of the classical trajecto- 
ries x(t) of complex Hamiltonians were taken up, 
and many interesting observations were reported. 
For example, it was found in some cases that the 
complex classical trajectories of a particle having 
real energy are closed and periodic, while those for 
a particle having complex energy are open and ir- 
regular 0, H, llj. However, this does not hold in 
general, as it was observed that for a special dis- 
crete set of curves in the complex-energy plane, the 
classical orbits are actually periodic, because of the 
possibility of some sort of quantization condition 
Q. Additionally, attempts have even been made 



to explain the concept of tunnelling with the help 
of classical trajectories |7|. 

However, all the previous studies of complex 
classical systems have been for bound states. In 
one of our previous studies, we had obtained the 
exact analytical classical trajectories of the bound 
states of the VT symmetric Scarf II potential @ . 
The main observation in rcf. [9[ was that the 
switching of energy values from real to complex 
ones at the VT transition point in the quantum 
picture, was accompanied by the changing of the 
corresponding classical trajectories from closed pe- 
riodic curves to open irregular ones. Naturally, 
our next attempt is to obtain the exact analytical 
classical trajectories of the scattering states of the 
same potential. In particular, our main aim in the 
present work is to see if spectral singularity in the 
quantum picture reveals any interesting feature in 
the corresponding classical version. For this pur- 
pose, we shall follow the factorization technique 
developed by Kuru et. al. [2l|. The motivation 
for studying the Scarf II potential is varied — 

1. The problem is exactly solvable both quantum 
mechanically as well as classically @. 

2. The bound state spectrum exhibits the inter- 
esting phenomenon of VT phase transition at the 
exceptional point, where the real spectrum enters 
the complex domain [22}. 

3. The scattering state spectrum shows a spectral 
singularity at a single positive energy, where both 
reflection and transmission coefficients diverge and 
the eigen states are no longer linearly independent 

[II, i|. 

4. This potential has been used to study nonlin- 
ear optical beam dynamics in a single VT complex 
crystal [14 1. 



The organization of the paper is as follows. In 
Section II, we briefly introduce the formalism ap- 
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plied in ref. [21], and extend the same to the 
complex domain. We apply it to obtain the ex- 
act analytical classical trajectories for the classical 
analogue of the real (Hermitian) Scarf II potential 
in Section III, and for the VT symmetric poten- 
tial in Section IV. Section V is specially devoted 
to the classical picture at the spectral singularity. 
Finally, Section VI is kept for Conclusions and Dis- 
cussions. 



so that A ± and H are assumed to define a de- 
formed algebra 

{A ± ,H} = ±ia{H)A ± , {A + , A~} = -i/3{H)A ± 

The unknown functions are determined by solving 
the following equations simultaneously 



a(H) 



<p'(x)+g'(x)VH 



(7) 



II. FORMALISM 

This section is primarily included to make the 
work self-contained. To briefly discuss the for- 
malism applied in ref. [2l|, we start with the 
one dimensional classical Hamiltonian (in units 
H=2m=l) 



H{x, P ) = p 2 + V{x) 



(1) 



where V(x) denotes the potential, and x and p 
are the canonical coordinates of position and mo- 
mentum, respectively, with their Poisson bracket 
{x,p} = 1. Hence, the classical particle obeys the 
equations of motions as given by the Hamilton's 
equations 



. 9H „ 

x = — = 2p 
op 



P = —^ = -V'(x) (2) 



so that x — 2p = 2V'(x) which on integration 
gives the velocity of the particle as 



dx 



V(x) 



(3) 



E being its energy. While extending the formalism 
to the complex domain, the particle is expected to 
lie in the complex plain, so that the path x(t) it 
traces out, as well as its velocity v(t), may take 
complex values. The initial conditions determine 
the initial velocity of the particle, and any point 
in the complexplane may be taken as an initial 
starting point @. 

Following the factorization technique of ref. 
plj , we assume a factorization of the Hamiltonian 
H in the form 



H = A+A- +j(H) 



(4) 



In usual quantum mechanical factorizations, 7 is 
the factorization constant, but in eq. (U]) above, 
j(H) may depend on H. Furthermore, A^ are no 
longer complex conjugates. are taken to be of 
the form 

A ± = T if( x )p + V~Hg(x) + <p(x) +4>{H) (5) 



f(x)V'(x)-2f'(x)[H-V(x)} 

a(H) {g{x)VH + ip{x) + (j>{H)] !S j 



(9) 



2VH {f(x)g(x) - f(x)g'(x)} 
-^=g{x) [2f'{x)V{x)+f{x)V'{x)] 



H 



+ lf{x)<i>'(H)[H-V{x)] 
2f{x)&{x)+<P'{H)V'(x)] 

Now we construct two time dependent integrals of 
motion of the form 

Q ± = A ± e^ la ^ t (10) 

so that the values of may be denoted by 

q ± = c{E)e ±l8 ° (11) 

where 9q is determined from initial conditions, and 

c(E) =| q ± |= y/E—Y(H) (12) 

For c(E) to be real, the expression within the 
square root sign must be positive. This condition 
gives the range of energy values for the classical 
particle. It is worth noting here that for bound 
and scattering states, the complex character of the 
factors A^ may change and the deformed algebra 
may also be different [2l| . 



III. CLASSICAL ANALOGUE OF THE 
REAL SCARF II POTENTIAL 

This section is included to provide a direct com- 
parison between the classical pictures of the corre- 
sponding Hermitian and VT symmetric Scarf II 
potentials. For bound states (i.e., E < 0), A ± are 
defined as 

A ± = Tif{x)p + \f^Hg{x) + V {x) + (t){-H) (13) 

For unbounded motion (E > 0), putting \J — E = 
i\J\ E |, changes the algebra of the Poisson bracket 
to 

{A ± ,A }=±i^A ± , {A + ,A-} =za A (14) 
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where Aq = —ivH. 

The final expression for V(x) should be of the form 

V (x) — 70 sech — - — h 26 sech — — tanh — — 

2 2 (15) 

We have intentionally considered the coefficient of 

2 OLqX . , . 

sech — ~ — . viz., 7q, to be positive, since the quan- 



tum mechanical VT symmetric version of eq. (jf 51) 
displays the interesting phenomenon of spectral 
singularity. We shall return to this point later, 
in Sections IV and V. Eq. (fl"5"j) . demands that in 
the expression for A ± in (TT31 we take the following 
forms of the functions g(x) , (p(x) , <j){H) : 



g(x) ^ , tp(x) = , </>{H) = -i 



H 



(16) 



Solving equations (0, dHJ and ((9|) simultaneously, 
one of the possible options could be 



g(x) = sinh ^ , f(x) = cosh ^ (17) 



Pit) 



cosh 



CtQX 



z{E) cosh (Oq + <XQ\fEt\ (25) 



where 6$ denotes the initial starting point. 
In Fig. 1, we plot the Hcrmitian Scarf II poten- 
tial as in eq. (|T5"j). and in Fig. 2, we plot the 
phase space trajectories for this potential, for dif- 
ferent positive energies. Parameter values used are 
ao = 2, (5 = 2, 70 = 6. The plot shows that as the 
particle comes near the barrier, its classical mo- 
mentum, and hence kinetic energy, decreases, as 
expected. 
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FIG. 1: Plot of V(x) against x for Real Scarf II pot. 



Putting 



a(H) = ia V~H 



l{H) = 7o + 



H 



(18) 



(19) 



in eq. (j4]), we obtain after some simplification, the 
Hermitian Scarf II potential 



2 a x 



aox 



(XqX 



H = p + 70 sech — - — h 26 sech—— tanh — — 



Hence 



2 
(20) 



h lor real Scarf II poUrntial. 1'or 




FIG. 2: Plot of p(t) against x(t) for classical analogue 
of Real Scarf II potential 
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c {E) = \jE- 10 -- 



(21) 



Since we are dealing with scattering states, E is 
always positive; therefore 



E > 



y/H + 4r5 2 + 7 o 



The integrals of motion Q ± , are pure imaginary; 
so values of take the form 



q ± = Ti c(E) exp[=F6> ] 



(23) 



Straightforward calculations give the classical tra- 
jectories and momenta for the classical analogue of 
the quantum mechanical Hermitian Scarf II model 
as 

x(t) = — sinh -1 sinh (<9 + a VEt] +41 

a { ^JE V / h J 

(24) 



IV. CLASSICAL ANALOGUE OF THE VT 
SYMMETRIC SCARF II POTENTIAL 

For VT symmetric Scarf II potential, we put 
5 = iSj, so that 



(22) H = p z + 7q sech 2 ^^ + i 2 <5/ sech^^ tanh 

2 (26) 



aox 



aox 



For scattering states, this particular quantum 
mechanical model exhibits the interesting phe- 
nomenon of spectral singularity at a single positive 
energy E s 



-(70 + ~ 



(27) 



Spectral singularity (ss) is also known as zero- 
width resonance, as both the transmission and re- 
flection coefficients diverge at this particular en- 
ergy. Thus these (ss) may be seen as positive en- 
ergy discrete poles of transmission and reflection 
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coefficients. For the potential given above in eq. 
(|26p . in the quantum picture, ss occurs when the 
parameters 70 and Sj satisfy the following condi- 
tions m mi 

I 2(5/ I > 70 + ^ 

3 

70+ I 2(5/ I = 4n 2 + An + - , n = 0, 1, 2, • • • 

The main motivation of the present study is to see 
if this feature is manifested in the corresponding 
classical picture in any way. 

For scattering states and hence positive energy 
with 



c(E) = ^E- la + ^ (29) 
the range of energy values is obtained as 

either < E < 70 = ^ = g* (30) 



p(t) = Re p(t) + i Im p(t) , etc. In Fig. 3, we plot 
the potential profile of the VT symmetric Scarf II 
Hamiltonian, for parameter values ao — 2, 5i = 
2, 70 = 6. In Figures 4, 5, 6 and 7, we plot a se- 
ries of classical trajectories and momenta for this 
particular model, for both real and complex en- 
ergy, taking the energy to lie within the permissible 
range, as given in equations (|5D|) and ([3"T]) . 



Phase space trajectories for PT sym Scarf II potential, for real, positive P 




Re Kit) 



FIG. 4: Plot of Re p(t) against Re x(t) for classical 
analogue of VT symmetric Scarf II potential, for real 
energy, when 70 > | 2(5/ | 



or E > WlHgl (31) 

If the real part of the potential is greater than the 
imaginary part, i.e., 70 > | 2(5/ |, then energy is 
always real. But in the reverse case, i.e., when the 
imaginary part of the potential is greater the real 
part (70 < I 2(5/ |) then energy becomes complex, 
with a positive real part. The interesting point to 
note here is that even beyond the VT threshold 
(or phase transition point), nonlinear states can 
still be found in the quantum version of the Scarf 
II potential, with real eigenvalues [14j . 



Fig. 4 shows the real part of the phase space trajec- 
tories for the classical analogue of the VT symmet- 
ric Scarf II potential, for same parameter values as 
in Fig. 3 : ao = 2, (5/ = 2, 70 = 6. Analogous to 
the Hermitian case shown in Fig. 2, the particle 
slows down at the potential barrier. In both Fig. 2 
and Fig. 4, the phase space curves are for different 
energies. 



Classical trajectories for PT symmetric Scarf I! potential, lor real, positive energy 
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FIG. 3: Plot of Re V(x) and Im V(x) against x for 
the VT symmetric Scarf II potential 

Proceeding in a straightforward way, we obtain 
expressions for the classical trajectories and mo- 
menta, similar to those in equations (|24"|) and (|25|) , 
with S = iSi , c(E) = Re c(E) + i Im c(E) 
E = E R + iEi , x{t) = Re x(t) + i Im x(t) 



FIG. 5: Plot of Im x(t) against Re x(t) for classical 
analogue of VT sym. Scarf II potential, for real energy 



Figures 5 and 6 show the classical trajectories 
traced out by the particle, for real and complex 
energy respectively. The parameter values are 
a = 2, (5/ = 2, 70 = 6 ( 7o > | 2(5/ |), for the 
former, and ao — 2, (5/ = 2, 70 = 3 (70 < | 2(5/ |) 
for the latter. The curves in Fig. 6 show a defi- 
nite distortion from those in Fig. 5. Evidentally, 
figures 5 and 6 have no counterpart in the classical 
analogue of the Hermitian Scarf II model. Similar 
to the bound state curves of ref. [9j], the initial 
starting point determines the trajectory, for the 
same energy. 
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Classical trajectories fur FT symmetric Scarf II potential, tor complex energy 




FIG. 6: Plot of Im x{t) against Re x(t) for classical 
analogue of VT sym. Scarf II potential, for complex 
energy 



V. CLASSICAL PICTURE AT SPECTRAL 
SINGULARITY 

It may be mentioned that spectral singularities 
(ss) are typical for complex scattering potentials, 
with no analogue in Hermitian systems. A spectral 
singularity of H occurs at a positive energy E s = 
k% in the continuous spectrum of H where the Jost 
solutions ipk±(x), satisfying ipk±( x ) e ±tkx for 
x — > ±00 , E — fc 2 , become linearly dependent 
and their Wronskian goes to zero [24j. Thus, phys- 
ically, spectral singularities correspond to scatter- 
ing states that behave like zero- width resonances. 
Spectral singularities may be seen as positive en- 
ergy discrete poles of transmission and reflection 
coefficients. Physical realization of spectral sin- 
gularities are possible in VT symmetric optical 
waveguides. If the frequency u of the incoming 
wave in such a waveguide can be tuned to the fre- 
quency uj s of the spectral singularity, then the am- 
plitude of the outgoing wave will be considerably 
enhanced as lo — > lj s . 

To study the phenomenon of spectral singular- 
ity in the classical framework, we plot the real part 
of the phase space trajectories for the VT sym- 
metric Scarf II potential, in Fig. 7, for different 
values of real energy, but when the real part of the 
potential is smaller than the imaginary part, i.e., 
70 < I 25 1 |. To be more precise, the plots are 
for parameter values 5i = 12, 70 =4, ao = 2, but 
for different positive energies. For this particular 
set of parameter values, it is observed that for low 
energies (as shown in the inset figures on below 
left), the classical particle is unable to overcome 
the barrier, and the phase space trajectories show 
a discontinuity near the potential barrier. As the 
energy increases the two disjoint curves come closer 
together, showing the inaccessible region for the 
classical particle decreases. Suddenly, as the en- 
ergy reaches a specific value viz., E s (13.7 for this 



particular set of parameter values), the real part 
of the classical momentum at the potential bar- 
rier tends to diverge. Beyond this energy, the mo- 
mentum decreases with increase in energy, finally 
reaching a constant value. As compared to the 
phase space trajectories in Fig. 4 (70 > | 25 1 |), 
these plots are inverted — the classical momen- 
tum decreases as the particle moves away from the 
potential barrier. We can identify this specific en- 
ergy E s as the energy at which classical spectral 
singularity takes place. 

Thus the anomalous behaviour of a particle at 
the spectral singularity is found in both the quan- 
tum and classical pictures, though the manifesta- 
tion is different in the two scenarios — while the 
reflection and transmission amplitudes blow up in 
the quantum picture, the real part of the momen- 
tum tends to diverge in the corresponding classi- 
cal picture. This is the most important finding of 
the present study. Another point worth noticing 
here is that the condition for spectral singularity 
is slightly different in the quantum and classical 
pictures. While the quantum condition for ss is 

I 25 1 I > 70 H — , the classical condition for the 

same turns out to be | 25 / | > 70. This finding 
is similar to the condition for exceptional point in 
the two pictures [9(. 



VI. CONCLUSIONS 

To conclude, we have studied the scattering 
states (having positive energy) of the classical ana- 
logue of an exactly solvable, non Hermitian quan- 
tum mechanical Hamiltonian, viz., the VT sym- 
metric Scarf II potential 

V(x) — 70 sech 2 a; — 2i5j sech x tanhx 

by applying the factorization technique of ref . [2l[ . 
In particular, we have obtained the exact analyt- 
ical classical trajectories of a particle subject to 
this potential, and moving about in the complex 
x plane. The main motivation behind this study 
was to see if the interesting phenomenon of spec- 
tral singularity exhibited in the quantum version 
at a particular energy (say E s ), is manifested in 
the classical picture as well. 

We plotted a series of curves, for different pa- 
rameter values, for a clear picture of our findings. 
Figures 1 and 3 show the Hermitian and VT sym- 
metric Scarf II potentials, respectively. Fig. 2 
gives the classical phase space trajectories for the 
former, while Fig. 4 gives the real part of the same 
for the latter. That the two figures resemble each 
other is proof of the fact that the complex VT 
symmetric potential shares the classical features 
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Phase space curves for PT sym Scarf II potn, at spectral singularity 
£,= 13.7 
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Colour Online : 



| Dashed (Green) curve for E > Es, 
j DotDashed (Black) curve for E » Es, 

Dotted (Blue) curve for E = Es, 
■ Solid (Red) curve for £ < Es, 
Magenta curve for E « Es, 
Cyan (bottom left) 
curve for least E « Es, 

Inset plots give better picture 




FIG. 7: Colour Online : Plot of Re p(t) against Re x(t) 
for classical analogue of VT symmetric Scarf II poten- 
tial, exhibiting spectral singularity at a single real en- 
ergy E a = 13.7, when 70 < | 2Si | ; Inset curves are 
shown for a better understanding of this phenomenon 
— the 3 plots on below left are for E < E s with the 
bottom most curve of least E, while the 3 on upper 
right are for E > E B with the topmost curve of highest 



of the real (Hermitian) potential. Fig. 5 and Fig. 
6 show the classical trajectories traced out by the 
particle (Im x{t) against Re x(t)), in the VT sym- 
metric potential, for real and complex energy, re- 
spectively. The point to be noted here is that none 
of the trajectories cross each other in any plot. 

Fig. 7 shows in detail the strange behaviour 
of the classical particle at the spectral singularity, 
when the imaginary part of the potential is greater 
than the real part, but energy is real. So long as 
the energy of the classical particle is sufficiently 
low, in fact lower than E s , it is unable to cross 
the barrier, and the phase space trajectories show 
a discontinuity near the potential barrier. All of 
a sudden, at the specific energy for spectral singu- 
larity E s , the real part of the classical momentum 
tends to diverge as the particle reaches the barrier. 
Thereafter, the particle momentum at the barrier 
decreases with increasing energy, thus showing a 
peculiar behaviour. 

To give more stress to the strange behaviour of 
the classical momentum near and at the spectral 
singularity (ss), and to see how one obtains the 
specific energy (E s ) at ss, we plot the real part of 
the classical momentum against energy, near the 
centre of the potential barrier, in Fig. 8, for the 
same set of parameter values as in Fig. 7 : Si = 
12, 70 = 4, ao = 2. For this particular set of 
parameter values, the spectral singularity is found 
to occur in the classical picture at E s — 13.7 . 
As the energy approaches E s , the momentum rises 
sharply, then falls rapidly with increasing energy, 
finally attaining a constant value for large energy. 

Plot of classical momentum vs energy 




Energy 



FIG. 8: Plot of Re p against Energy for classical ana- 
logue of VT symmetric Scarf II potential, exhibiting 
spectral singularity at a single real energy E a = 13.7, 
when 81 — 12, 70 = 4, «o = 2, i.e., 70 < | 25/ | 



In our earlier work on bound states of the VT 
symmetric Scarf II potential 9] , we had observed 
that at the VT transition point, also known as 
the exceptional point, when the energy spectrum 



E 
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in the quantum picture goes from the real to the 
complex domain, the classical trajectories switch 
from closed, periodic orbits to open, irregular ones. 
In the present study on the scattering states for 
the same potential, we observe that at the spec- 
tral singularity, when reflection and transmission 
coefficients blow up in the quantum version, the 
classical momentum tends to diverge. Addition- 
ally, similar to the condition for exceptional point 
0, the classical condition for spectral singularity 
differs from the corresponding quantum condition 

by a factor of -. 
J 4 
Our next attempt would be to study more such 

cases, to check whether this strange behaviour of 



the classical particle at the spectral singularity is 
characteristic of this particular potential, or some 
generalized conclusion can be drawn regarding the 
classical behaviour at the spectral singularity. 
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